Josephson oscillations between exciton condensates in electrostatic traps 
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Technological advances allow for tunable lateral confinement of cold dipolar excitons in coupled 
quantum wells. We consider theoretically the Josephson effect between exciton condensates in two 
traps separated by a weak link. The flow of the exciton supercurrent is driven by the dipole energy 
difference between the traps. The Josephson oscillations may be observed after ensemble average 
of the time correlation of photons separately emitted from the two traps. The fringe visibility is 
controlled by the trap coupling and is robust against quantum and thermal fluctuations. 
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I. INTRODUCTION 

The Josephson effect is a a macroscopic coherent phe- 
nomenon which has been observed in systems as diverse 
as superconductors^ superfluid Helium; 2 . Bose-Einstein 
condensates in trapped ultra-cold atomic gases 4 Since 
Josephson oscillations appear naturally when two spa- 
tially separated macroscopic wave functions are weakly 
coupled, they have been predicted for bosonic excitations 
in solids as well, like polaritons^ and excitons.- How- 
ever, unlike the polaritons, which have a photonic com- 
ponent allowing for easy detection^ excitons stay dark 
unless they recombine radiatively. So far, it is unclear 
how the exciton Josephson effect could be observed. In 
this paper, we propose an experiment. 

Condensed excitons are predicted to emit coherent 
light. 8 If Josephson oscillations occur between two exci- 
ton traps, in principle they can be probed by measuring 
the interference of the beams separately emitted from the 
traps. However, in the time interval before recombina- 
tion, there are too few photons emitted for an adequate 
signal to noise ratio, and one has to average the signal 
over many replicas of the same experiment. 9 . We will show 
that such ensemble averaging blurs the signature of the 
Josephson effect except in the relevant case of exciton 
"plasma" oscillations. 10 For the latter the dipole energy 
difference between the traps modulates the visibility a of 
interference fringes, providing a means for detection. 

The paper is organized as follows: In Sec. |TT] we in- 
troduce the double quantum well system and illustrate 
a feasible scheme to manipulate electrically the exciton 
phase. After setting the theoretical framework (Sec. lIIIjl . 
we discuss the proposed correlated photon counting ex- 
periment (Sec. IIVP and provide an estimate for its key 
parameters (Sec. fVT). 



II. ELECTRICAL CONTROL OF THE EXCITON 
PHASE 

Consider a double quantum well where electrons and 
holes are separately confined in the two layers. In ex- 
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FIG. 1: (a) Double quantum well energy profile along the 
growth direction z. (b) Exciton potential profile of the double- 
trap system along x. (c) Arrangement for measuring the time 
correlation of the emitted photons from the two traps. 



periments aiming at Bose-Einstein condensation of ex- 
citons, electron-hole pairs are optically generated off- 
resonance, left to thermalize, form excitons, and, at suf- 
ficiently low temperature T and high density, condense 
before radiative decay.— Let z be the growth axis of 
the two wells separated by distance d. The electrons 
in the conduction band and holes in the valence band 
move in the planes z = d, 0, respectively [Fig. [If a)]. Let 
^ a (x, y, 0, t) and \&£(x, y, d, t) denote the respective elec- 
tron and hole destruction operators, with the vacuum 
being the semiconductor ground state with no excitons. 
In the experiments^ an electric field F z is applied along 
z to suppress inter-layer tunneling, thereby quenching 
the exciton recombination. Fabrications^ 2 - of electrostatic 
traps with suitably located electrodes to provide lateral 
confinement for the excitons have been implemented. 
The double quantum well is sandwiched between two 
spacer layers, providing insulation from planar electrodes 
lithographed on both sides of the coupled structure. Each 
electrode controls a tunable gate voltage, V g (x,y,z), 
which localizes in a region of the xy plane the field com- 
ponent along z, F z (x,y,z) — —dV g (x,y,z)/dz, while F x 
and F y are small and can be neglected as well as the the 
dependence of F z on z. The vertical field F z (x,y) makes 
the electrostatic potential energy of the exciton dipole 
depend on the lateral position, Ux(x,y) — —edF z (x,y) 
(e < 0) [cf. Fig. Hfb)]. In this way, potential traps for 
excitons are designed with great flexibily, with in situ 
control of the height, width, and shape of the potential 
barriers^ 

First, we focus on the quasi-equilibrium situation be- 
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fore radiative recombination, where excitons condense in 
two coupled electrostatic traps, both within the conden- 
sate coherence length. Figure QJb) depicts the exciton 
potential profile Ux{x,y — 0) along the x axis, with a 
link between two identical traps. The potential barrier 
allows tunneling between the condensates Si (x, y, t) and 
S2(x, y, t) in the two trapsJ^ The optical coherence in a 
single trap is of the form 



S(a, y, t) = (*t(a:, y, 0, t)* 6 (ar, y, d, t)) 



(1) 



where (. . .) denotes quantum and thermal average. In 
the limit nog <C 1, with ob being the two-dimensional 
effective Bohr radius and n the exciton density, 3(ar, y, t) 
is the macroscopic wave function for the center-of-mass 
motion of excitonSf« which may be written in the form 



2{x,y,t) = y/n 



(2) 



with n s being the density of the exciton condensate and 
ip the phased For a gauge transformation of the gate 
potential V g — ► V g — c~ 1 dx{t)/dt, which leaves the field 
F z unaltered, the field operators iff gain a phase, 



repulsive character of the dipolar interaction between ex- 
citons in coupled quantum wells. The condensate total 
wave function solution is 



£(x, y, t) = E 1 (x, y, Nje* 1 + E 2 (x, y, N 2 ) e** 



(6) 



where both the trap population Ni(t) and the condensate 
phase <pi(t) possess the entire time dependence for the ith 
trap (i = 1, 2), and 3$ (a:, y, Ni) is a real quantity, with 



dx dyZ*(x,y ) N i )=N i (t). 



(7) 



The dynamics of the GP macroscopic wave function 
3(x, y, t) depends entirely on the temporal evolution of 
two variables, the population imbalance k(t) — (Ni — 
N 2 )/2 and the relative phase <j>(t) = ipi — ip 2 of the 
two condensates. Here we consider a time interval much 
shorter than the exciton lifetime (10 — 100 ns) and ig- 
nore the spin structure. Therefore, the total population 
is approximately constant, Ni(t)+N 2 (t) = N. The equa- 
tions of motion for the canonically coniugated variables 
hk and d> are derived from the effective Hamiltonian 



iff a -> * Q exp 
* h -> * fe exp 
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x( x ,y,d, t) 



(3) 



The macroscopic wave function, by Eq. |T]), also gains a 
phase, 



V + ir [x{z = d,t)-x(z = 0,t)] 
he 



(4) 



Hence, the frequency of time oscillation of the condensate 
is given by the electrostatic energy of the exciton dipole 
in the external field, 16 U — —edF z : 



(p = tp 



1 
h 



(0) _i_ ^ edFzt 



(5) 



with ip(°^ being the time-independent zero-field valued 
In the absence of the bilayer separation of the electrons 
and the holes, their gauge phases gained in the electric 
field would cancel each other resulting in no time depen- 
dence driven by U. Equation ([5]) shows that the exper- 
imentally controllable dipole energy difference between 
the two traps depicted in Fig. QJb), AU = —ed(F z i — 
F z2 ), drives the relative phase between the two conden- 
sates, thereby creating Josephson oscillations as a means 
for measuring the Josephson tunnel between the traps. 



III. EXCITON JOSEPHSON OSCILLATIONS 

We next introduce the usual two-mode description 
of inter-trap dynamics based on the Gross-Pitaevskii 
(GP) equationj^££ i 10 i 18 i 19 Exciton-exciton correlation^ 
beyond the GP mean field may be neglected due to the 



Hj 



k 2 



AUk - -^-\/N 2 -Ak 2 cose 



(8) 



under the condition k < N (Ref. 10). E c = 2dfi\/dN\ is 
the exciton "charging" energy of one trap, where [i\ is the 
chemical potential of trap 1, whereas Sj is the Bardeen 
single-particle tunnelling energy, 



Sj = — / dy 

m 



(9) 



x=0 



where m is the exciton mass. The single-particle orbital 
£i(x,y) is defined through Ei(x,y) = y/Ni£i(x,y). 

The various dynamical regimes associated to certain 
intitial conditions (k(0),<p(0)), including tt oscillations 
and macroscopic quan tum self-trapping, are exhaustively 
discussed in Refs. \fQ Two cases are specially relevant: 



1. AC Josephson effect 

Under the conditions AU > NE c /2, AU > 5j, one 
easily obtains 



k = sin ( 

2h 



(10) 



Equation (|10p shows that, analogously to the case of two 
superconductors separated by a thin barrier, if the phase 
difference 4> between the condensates is not a multiple of 
7r, an exciton supercurrent 2k flows across the barrier. 
Remarkably, in the presence of an electric field gradient 
along z, an exciton flux oscillates back and forth between 
the two traps, with frequency AU/K. As an exciton goes 
through the barrier, it exchanges with the field the dipole 
energy acquired or lost in the tunneling process. The 
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analogy with the AC Josephson effect for superconduc- 
tors is clear: in that case a bias voltage V is applied 
across the junction, and the energy 2eV is exchanged be- 
tween field and Cooper pairs, as the latter experience a 
potential difference of V when penetrating the potential 
barrier. 



2. Plasma oscillations 

This case concerns small oscillations around the equi- 
librium position (fc, (/)) eq = (0,0). The Hamiltonian ([8]) 
may then be linearized into the form 

It follows that both k and <j) oscillate in time with plasma 
frequency 



ojj = -y/5j(NE c /2 + Sj). (12) 

Note that AU displaces the equilibrium position from 
(M)eq = (0,0) to 

(k, (£) eq = (-AUN5j/2(hLuj) 2 , 0). (13) 

IV. CORRELATED PHOTON COUNTING 
EXPERIMENT 

Figure [TJc) illustrates the correlated photon count- 
ing setup which we propose to probe Josephson oscil- 
lations. The detector measures the intensity I(t) of 
the sum of the two beams separately emitted from the 
traps. A delay time r is induced in one of the two 
beams, as in Ref. § The fields are simply propor- 
tional to the order parameters E< of the traps. In fact, 
H(x, y, t) is associated with a macroscopic electric dipole 
moment, Pit) = xP x (t) ± iyP y (t), which couples to 
photons: P x (t) — J dx dy x S(x, y, t), and similarly for 
P y . The built-in dipole (Pit)) ^ oscillates with fre- 
quency ifi + Ex)/h, where Ex is the optical gap minus 
the exciton binding energy, and fi accounts for exciton- 
exciton interaction.— This macroscopic oscillating dipole 
is equivalent to a noiseless current, which radiates a co- 
herent fieldJi Therefore, the measured intensity J(r) is 
7(r) = 2I [1 + (cos</)(r))], assuming that the fields emit- 
ted from the two traps have the same magnitude (and 
intensity 1$) but different relative phase (f>, which is eval- 
uated at the delayed time I(r) may be written as 

/(r) =2/ o [l + acos0o(r)] , (14) 

where (/)o( T ) is the phase averaged over many measure- 
ments, defined by the condition (sin [<p(r) — </>o(t)]) = 0, 
and 

a = (cos[0(r)-0 o (r)]) (15) 



is the fringe visibility i.e., the normalized peak-to- valley 
ratio of fringes, a = (7 roax -imin)/ (4ax+^min), with 7 max 
(Jmin) being the maximum (minimum) value of 7(r), and 

< a < 1. 

Equation (TT4"]) has a few important caveats. Since /(r) 
is an average, the temporal inhomogeneous effect will 
blur the interference fringes, i.e., a < 1. Other dephas- 
ing mechanisms include exciton recombination and in- 
elastic exciton-phonon scattering^ as well as inelastic 9 
and elastic*^ exciton-exciton scattering, which in first 
instance may all be neglected for short t, low T, and 
nag <C 1 , respectively. The most immediate caveat is 
that the exciton condensates in decoupled traps must 
acquire a relative phase if initially they condense sep- 
arately without a definite phase relation. This scenario 
is analogous to the case of interference between indepen- 
dent laser sources first discussed by Glauber— and later 
studied experimentally for matter waves; 2 ^ Even though 
a one-shot measurement with sufficient resolution would 
display fringes, the relative phase </>o( r ) is also subject 
to intrinsic dephasing effects by quantum fluctuations. 21 
The latter are significant noise sources which affect a, 
when <f> and k are quantized into canonically conjugated 
quantum variables whereas in the CP theory used so far 
they were classical variables whose fluctuations where ne- 
glected. 

In the following, we quantize Hamiltonian ([8]) in or- 
der to properly evaluate a = (cos [0(r) — <^o(r)]) as a 
quantum statistical average in finite traps. Therefore, we 
follow Ref. and introduce the commutator [4>, k] = i. 
The operator k now appearing in the quantized version 
of Hamiltonian ([H]) takes the form —id/d(f>, whereas the 
ground state wave function is defined in the space of pe- 
riodical functions of <fi with period 2ir. If condensate os- 
cillations are mainly coherent, the variance of </> is small 
and the visibility is approximated by a — 1 — h ((A</>) 2 ^. 

The most interesting case concerns plasma oscillations. 
For AU — 0, the ground state of the quantized version 
of the harmonic oscillator Hamiltonian (jTTj) is a Gaus- 
sian, with 4>q = 0, independent from r, and minimal 
spreading (Aq9) w iE c /25. 1 N) 1 / 2 . Therefore, the inter- 
ferometer output is time- independent, / = 2/o(l + ®)-, 
showing constructive interference, I > 21$, with a = 

1 — iEcl&SjN) 1 / 2 . Not surprisingly the visibility is 
controlled by the ratio E c /5jN, reaching the maximum 
a = 1 as E c /SjN — > 0. In fact, a is given by the balance 
between the competing effects of tunnelling (cx 5jN), 
which enforces a well-defined inter-trap phase, and in- 
verse compressibility (cx E c ), which favors the formation 
of separate number states in the two traps, thus separat- 
ing the two macroscopic wave functions. 

A small finite value of AU in Eq. (TTTj) displaces the 
equilibrium position of the harmonic oscillator. Notice- 
ably, the ground state is a coherent state with non-null 
evolution of the average phase in time, 

<Mt) = - sin (ujt), (16) 
Tito j 



FIG. 2: (Color online). Beam intensity I{t)/Iq vs delay time 
r, for AU/huij = 0.2, 0.5 (dashed and solid lines, respec- 
tively), (a) T = and a = 1, 0.8 (black and red [light gray] 
lines, respectively), (b) a(T = 0) = 0.94 and fesT/Euj = 
0, 1, 2 (black, red [light gray], and blue [dark gray] lines, 
respectively) . 



whereas a is unchanged. This key feature allows 
for directly monitoring r-dependent plasma oscillations 
of frequency ljj through the photon correlation mea- 
surement (cf. Fig. [2]). We evaluate the effect of 
thermal fluctuations on a via the formula a(T) — 
Y, n (Xnexp[-l3E n ]/Y, n fxp[-l3E n ], where (i = l/k B T, 
ks is the Boltzmann constant, 2(a n — 1) = ((A0) 2 ) n is 
the variance of </> in the nth excited state whose energy is 
E n . At low T, the excited states may be approximated 
as those of the harmonic oscillator, giving 



a{T) = 1 



E r . 



1 



28 jN \2 ' e^-> - 1 
The above results are summarized by the formula 



(17) 



I(r) = 2J 



1 + a(T) cosf — — sinwj r 



(18) 



which is valid for E c /6jN <C 1. For small dipole 
energy variations, AU/fruij <C 1, the oscillating part 
within the square brackets of Eq. (| 18|) may be written 
as -a(T)/2(AU/hujj) 2 sin 2 ujt. This shows that the 
visibility a(T) of fringes, which oscillate like sin 2 uijt, 
is modulated by the experimentally tunable factor 
{AU/hujj) 2 /2. The dependence of 7(r) on AU is illus- 
trated in Fig.[2]for two values of AU /Tiloj. As AU /Tiuij is 
increased [from 0.2 (dashed lines) to 0.5 (solid lines)], the 
amplitude of oscillations of J(t) shows a strong non- linear 
enhancement, providing a clear signature of Josephson 
oscillations. The oscillation amplitudes are larger for 
higher values of a [cf. Fig.^a)], and fairly robust against 
thermal smearing [cf. Fig.[2{b)]. In fact, Fig.[^b) shows 
that the oscillation of 7(r) is still clearly resolved for 
temperatures as high as T k Tiujj/kB- At even higher 
temperatures a(T) displays anharmonic effects^ with 
a{T) — > as T — > oo. 



FIG. 3: (Color online). Beam intensity 7(r)//o vs delay time 
r, for AU/huj = 4, 12 (dashed and solid lines, respectively) 
and kBT/hw.j — 0, 2 (black [dark gray] and red [light gray] 
lines, respectively), with a(T = 0) = 0.94. 



The maximum peak-to-valley ratio of fringes attain- 
able for plasma oscillations is limited by the condition 
that AU shifts the equilibrium position (k,(f>) cq of the 
oscillator, as given in Eq. (fT5]) . slightly with respect to 
the origin: AU Sj/(hujj) 2 <C 1. For example, by tak- 
ing the values of ujj and Sj estimated in Sec. [V] and 
imposing the condition AU 5j/(hujj) 2 — 5 • 10 -2 , one 
has AU/Tiloj w 12. Figure [3] displays I{f) vs. r for 
AU /htuj = 4, 12 (dashed and solid lines, respectively). In 
both cases the range of amplitude oscillations of I is very 
close to the ideal interval [0, 4io] ■ Therefore, the intensity 
oscillations should be easily detected, even at finite tem- 
peratures (cf. the black [dark gray] and red [light gray] 
lines, corrsponding to ksT/fiDj = 0, 2, respectively). 
For large values of AU/hujj higher overtones appear in 
the oscillations of J(r) (solid lines in Fig. [3]), in addition 
to the fundamental frequency u>j, which is present for 
any finite value of AU (cf. dashed lines in Fig. [3]). 



V. ESTIMATE OF THE VISIBILITY AND 
PLASMA FREQUENCY 

We assess the feasibility of the experiment by estimat- 
ing the parameters of Eq. (fT8|) . Both a and ujj depend 
on E c and 5j. We evaluate the latter by first solving 
the GP equation for a two-dimensional harmonic trap 
within the Thomas-Fermi approximation^ and then by 
matching the wave functions of the two traps by using 
the semiclassical method of Ref. lltj. The coupling con- 
stant g — 4irde 2 /s r appearing in the non- linear term of 
the GP equation, multiplied by n s , is the energy shift of 
an exciton added to a parallel plate capacitor with sur- 
face charge density en s (s r is the quantum well dielectric 
constant) We obtain 



E c = 2huj Q N x 1/2 f — 

VaB 



1/2 



(19) 
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where Tiloq is the energy quantum of the trap and ag = 
h 2 e r /me 2 , as well as 

(20) 

Here u = 0.397, and S ~ 2 1 / 2 7r(V - fi)/fUD for a inter- 
trap quartic barrier of height Vo, with Vb — (J, <C Vb 
[Fig. [Tib)]. At density n = 2.5 • 10 10 cm" 2 , evaluated 
at the trap center, excitons are still weakly interacting 
(nag ~ 0.1 with ae ~ 20 nm). By taking GaAs parame- 
ters, d = 12 nm, N% = 10 3 , one has Uuoq — 11 /ieV, \i\ = 
440 fJ.eV, and a condensate radius of 1.6 fim. The barrier 
height Vb, as well as So, should be as low as possible. 
For Sq = 1 we obtain high visibility [a = 0.94 at T = 0, 
cf. Fig. [2][b)], as well as a plasma frequency lo,j/2tt of 
0.41 GHz, whose period (w 2 ns) is an order of magni- 
tude shorter than the exciton lifetime. Note that Tiloj — 
1.7 /xeV <C huio, hence the plasma oscillation is decoupled 
from single-trap modest The temperature associated to 
Tiljj, T = 20 mK, is very low but within experimental 
reach. 

The AC Josephson effect cannot be observed within 
our scheme. In fact, for large values of AU, the term 



proportional to cos</> appearing in the Hamiltonian (JSJ) 
may be neglected in first approximation, and the ground 
state wave function is a plane wave, (27r) -1 / 2 exp [in<j>], 
where n is the integer closest to —AU/E C . Since the 
probability density, (2tt)~ 1 , is constant, the phase is dis- 
tributed randomly and the visibility is zero. Therefore, 
the correction to a coming from the inclusion in the cal- 
culation of the term neglected in will be small and 
fragile against fluctuations. 



VI. CONCLUSION 

In conclusion, exciton plasma oscillations may be mea- 
sured by the time correlation of photon emission from two 
sides of the Josephson junction through electrical control 
of fringe visibility. Our findings pave the way to the ob- 
servation of the exciton Josephson effect. 
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